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Morphology  control  is  crucial  both  for  zinc— air  batteries  and  for  zinc— air  fuel  cells  during  zinc  regen¬ 
eration.  Zinc  dendrite  should  be  avoided  in  zinc— air  batteries  and  zinc  pellets  are  yearned  to  be  formed 
for  zinc— air  fuel  cells.  This  paper  is  mainly  to  analyze  the  mechanism  of  shape  change  and  to  control  the 
zinc  morphology  during  charge.  A  numerical  three-dimensional  model  for  zinc  regeneration  is  estab¬ 
lished  with  COMSOL  software  on  the  basis  of  ionic  transport  theory  and  electrode  reaction  electro¬ 
chemistry,  and  some  experiments  of  zinc  regeneration  are  carried  out.  The  deposition  process  is 
qualitatively  analyzed  by  the  kinetics  Monte  Carlo  method  to  study  the  morphological  change  from  the 
electrocrystallization  point  of  view.  Morphological  evolution  of  deposited  zinc  under  different  conditions 
of  direct  currents  and  pulse  currents  is  also  investigated  by  simulation.  The  simulation  shows  that 
parametric  variables  of  the  flowing  electrolyte,  the  surface  roughness  and  the  structure  of  the  electrode, 
the  charging  current  and  mode  affect  morphological  evolution.  The  uniform  morphology  of  deposited 
zinc  is  attained  at  low  current,  pulsating  current  or  hydrodynamic  electrolyte,  and  granular  morphology 
is  obtained  by  means  of  an  electrode  of  discrete  columnar  structure  in  combination  with  high  current 
and  flowing  electrolyte. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Appropriate  energy  storage  plays  a  more  important  role  in  a 
wide  range  of  applications,  including  frequency  regulation  of 
renewable  energy  sources,  peak-shaving  of  electrical  grid,  emer¬ 
gency  power  and  portable  power  1,2  .  Metal-air  batteries  have 
received  much  attention  due  to  their  high  energy  densities  [3], 
where  oxygen  from  the  air  employed  as  their  cathode  material. 
Compared  with  other  metals,  zinc  possesses  a  variety  of  advan¬ 
tages:  low  equilibrium  potential,  high  specific  energy,  good  envi¬ 
ronmental  compatibility,  good  reversibility,  it  is  found  in 
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abundance  and  it  is  easy  to  handle  [4] ;  therefore,  zinc  material  is 
widely  applied  to  the  electrochemistry  field  in  either  the  galvani¬ 
zation  or  in  cathode  active  material  5,6].  However,  several  major 
problems  in  zinc  regeneration  need  to  be  addressed  including  zinc 
pellets  of  appropriate  size,  shape  change  of  zinc  electrode  and 
dendritic  growth  [7].  These  are  particularly  important  for  zinc-air 
battery  due  to  seriously  affecting  the  battery's  durability  and  per¬ 
formance,  limiting  the  applications  of  the  zinc-air  secondary  bat¬ 
tery.  Moreover,  zinc-air  fuel  cells'  performance  is  subject  to  zinc 
pellets  size,  because  it  can  determine  active  surface  area,  and  in¬ 
fluence  electrolyte  flowing  and  contact  resistance. 

Although  such  morphological  changes  arise  from  ionic  trans¬ 
port  through  diffusion,  migration  and  convection,  the  mechanism 
for  causing  morphological  changes  remains  unknown.  Some  re¬ 
searchers  have  made  significant  progress  in  exploring  the 


66 


K.  Wang  et  al.  /  Journal  of  Power  Sources  271  (2014)  65—75 


mechanism  for  the  morphological  changes  of  zinc  electrodes 
[8-10].  McBreen  [11  developed  a  concentration-cell  model  that 
accounted  for  the  shape  change  of  zinc  electrodes  during  charge/ 
discharge,  but  for  slow  or  fast  diffusion  further  explanation  are 
needed.  Choi  et  al.  [12,13]  proposed  a  hydrodynamic  model  for 
shape  change  based  on  electrolyte  flows  in  a  battery  as  a  conse¬ 
quence  of  osmotic  and  electro-osmotic  forces  which  explains  the 
ionic  movement  from  the  zinc  electrode's  top  toward  its  bottom, 
but  this  method  cannot  provide  reasonable  explanations  for  the 
zinc  accumulation  at  the  center  of  the  electrode.  Bockris  et  al.  held 
that  dendritic  growth  might  be  closely  related  to  the  overpotential 
at  the  tip  of  the  dendrites  [8  .  They  also  investigated  the  rela¬ 
tionship  between  the  dendritic  formation  and  the  physical  prop¬ 
erties  of  an  electrolyte  such  as  composition,  temperature, 
viscosity,  and  hydrodynamics  [14  .  In  addition,  some  researchers 
have  done  some  work  in  solutions  to  zinc  dendrites.  Parker  et  al. 
[15]  designed  a  porous  electrode  to  suppress  growth  of  zinc 
dendrites.  Smedley  et  al.  16]  provided  a  solution  of  zinc  regen¬ 
eration  with  optimizing  the  structure  of  zinc  electrode  to  control 
zinc  morphology.  Shaigan  et  al.  17]  employed  pulsating  current  as 
morphological  control  of  electrodeposited  zinc.  Other  researchers 
have  also  made  some  attempt  to  model  the  zinc  electrode.  Deiss 
et  al.  18]  developed  a  one-dimensional  numerical  model  of  an 
electrically  rechargeable  alkaline  zinc-air  battery  that  can  predict 
decrease  of  the  zinc  redistribution  with  increasing  cycle  number. 
Song  et  al.  [19]  established  a  one-dimensional  mathematical 
model  of  the  discharge  process  and  of  the  failure  mechanisms  of 
zinc  electrode,  predicting  variable  distributions  within  the  elec¬ 
trode  itself.  Bauer  et  al.  [20 !  put  forward  a  model  of  the  numerical 
simulation  of  copper  electrodeposition.  Braff  et  al.  [21,22]  estab¬ 
lished  a  two-dimensional  model  for  simulating  a  hydrogen 
bromine  laminar  flow  fuel  cell  based  on  COMSOL  software.  To 
control  the  morphological  change,  the  mechanism  for  shape 
change  must  be  further  studied,  and  then  put  forward  reasonable 
solutions. 

In  this  work,  we  investigated  the  morphological  change  of  zinc 
electrode  related  to  parametric  conditions  with  COMSOL  software, 
and  experiments  on  zinc  electrodeposition  under  the  conditions  of 
different  current  and  flowing  electrode  were  also  conducted. 
Furthermore,  we  implemented  zinc  depositing  process  simulation 
relying  on  the  three-dimensional  kinetics  Monte  Carlo  method, 
analyzing  morphological  evolution. 

2.  Experiment  methods  and  numerical  model 

2.2.  Experiment  methods 


Table  1 

Testing  parameters  for  zinc  electrodeposition. 


Parameters 

Feature  I 

Feature  II 

Electrode  structure 

Planar  electrode 

Porous  electrode 

Electrolyte  state 

Quiescent 

Hydrodynamic 

Charging  mode 

Direct  current 

Pulsating  current 

2.2.  Numerical  model 
22.1.  Electrode  kinetics 

Zinc  electrodeposition  involves  energy  conversion  transforming 
electrical  energy  into  chemical  storage  energy,  accompanied  by 
phase  changes.  Fig.  1  shows  that  chemistry  and  electrochemistry 
reactions  occur  near  or  on  the  electrodes  of  a  zinc-air  rechargeable 
battery  during  charge,  where  the  morphological  change  of  zinc 
regeneration  is  a  concern  because  it  reduces  the  capacity  of  the 
battery.  Pulsating  currents  from  new  clean  energy  is  more  benefi¬ 
cial  to  suppress  dendritic  growth. 

Zinc  oxide  or  zincate  is  reduced  at  the  cathode,  and  hydroxyl  is 
oxidized  at  the  anode  during  charge.  The  electrochemical  reactions 
can  be  expressed  as: 

Anode:  20H~^io2T  +H20  +  2e  (1) 


Cathode  : 


ZnO  +  H20  +  20H"  -^Zn(OH)42“ 
Zn(OH)42~  +  2e  —  Zn  +  40H“ 


The  required  charge  voltage  can  be  determined  by  the  ther¬ 
modynamic  equilibrium  voltages  and  the  activation  overpotentials 
at  the  anode  and  cathode,  and  ohmic  loss  of  the  electrolyte.  Thus. 


A<p  —  Eg  q  —  Eg  q  +  Va  +  Vc  +  *^el  (3) 

Butler-Volmer  equations  are  used  to  describe  the  electrode 
reaction  kinetics. 


*o 


where  i\oc  denotes  the  local  current  density  at  the  solid-liquid 
interface,  i'o  exchange  current  density,  aa  anodic  transfer  coefficient, 
ac  cathodic  transfer  coefficient,  F/RT  the  Nernst  factor,  and  77  acti¬ 
vation  overpotential  for  the  electrode  reactions,  which  can  be 
defined  as  the  below  equation: 


An  experimental  apparatus  for  zinc  electrodeposition  was 
fabricated,  where  the  cathode  employed  stainless  steel,  the  anode 
used  nickel  network  with  the  same  dimension  as  the  cathode,  and 
the  electrolyte  was  made  of  weight  40%  potassium  hydroxyl  mixed 
with  zinc  oxide.  The  distance  between  the  anode  and  the  cathode  is 
5  mm,  and  the  effective  dimension  of  the  cathode  is  about 
30  mm  x  30  mm.  Zinc  electrodeposition  experiments  were  carried 
out  at  the  room  temperature,  where  the  power  supply  with  pul¬ 
sating  and  galvanostatic  mode  was  EA-PS  8065-10,  and  data 
acquisition  instrument  used  for  recording  current  signals  was  NI- 
PXI-1033.  In  addition,  auxiliary  components  were  not  described 
here  including  a  pump  used  for  the  flowing  electrolyte.  The 
morphology  of  deposited  zinc  was  examined  by  scanning  electron 
microscopy  (JSM-6460).  For  some  reason,  hydrodynamic  condi¬ 
tions,  charging  mode,  and  electrode  structure  can  influence 
morphological  change  of  electrodeposited  zinc.  Table  1  lists  the 
different  testing  parameters  including  flowing  electrolyte,  charging 
mode  and  structure. 


V  =  <Ps  -  </>/  -  £eq  (5) 

where  </>s  is  electrode  potential,  and  </>/  is  electrolyte  potential.  Ac¬ 
cording  to  the  Nernst  equation,  equilibrium  potentials  Eeq  are 
calculated  by: 


Anode 


Cathode 

Fig.  1.  Schematic  view  of  zinc  electrodeposition. 
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£e,(  Z,^/Zn),£»(Zn"/z„)-gl„^- 

Znz+ 


The  potential  drop  in  the  electrolyte  is  also  part  of  the  charging 
voltage,  which  depends  on  the  ionic  current,  electrolyte  concen¬ 
tration,  and  conductivity.  An  empirical  relation  can  be  described  by: 


equations,  four  field  variables,  and  nine  boundary  conditions.  We 
make  the  assumption  that  species  concentration  of  the  electrolyte 
is  well-distributed  at  the  flowing  direction  and  the  rest  parts  are 
insulated  in  addition  to  the  electrolyte  and  the  electrodes.  Table  2 
lists  the  relevant  variables  and  expressions  in  the  simulation  pro¬ 
cess.  Therein  the  electrolyte  conductivity  is  assumed  to  be  constant 
due  to  the  flowing  electrolyte  counteracting  concentration  change 
during  electrodeposition. 

3.  Results  and  discussion 


i  =  -crV^/  (7) 

where  a  is  the  electrolyte  conductivity,  and  </>/  is  the  potential  in  the 
electrolyte. 

2.2.2.  Species  transport 

Transport  of  ions  in  the  electrolyte  is  controlled  by  diffusion, 
migration,  and  convection.  Diffusion  stems  from  concentration 
difference,  migration  is  determined  by  additional  power,  and  con¬ 
vection  comes  from  fluid  dynamics.  The  material  balance  equation 
can  be  expressed  as: 


-7 jp  +  vVck  -  V-(DkVck  +  zkpkFckVcf)i )  =  Rks  (8) 

where  ck  stands  for  the  species  concentration,  v  the  fluid  velocity, 
Dk  the  diffusion  coefficient,  zk  the  charge  for  ionic  species,  fik  the 
mobility  of  the  charged  species,  subscript  k  designates  species  (i.e. 
Zn(OH)42^ ,  OFT ,  K+ ),  F  Faraday's  constant,  and  RktS  is  reaction  rate. 

The  model  of  electrolyte  movement  employs  incompressible 
Navier-Stokes  equation  as  follows. 

(dv  \ 

—  + wVvj  = -Vp  +  pV2v +f  (9) 

The  electroneutrality  condition  is  also  taken  into  account. 

ifzk-ck  =  0  (10) 

k=  1 

The  electrode  boundaries  will  move  during  zinc  deposition.  The 
molar  flux  of  zinc  Nzn  (mol  (m2  s)-1)  over  the  electrode  can  be 
described  by  Faraday's  law  according  to: 


^Zn  hoc 

nF 


(ii) 


where  the  stoichiometric  coefficients  for  zinc  is  vzn  =  1  for  both  the 
electrode  reactions,  and  the  number  of  electrons  is  n  =  2. 

The  depositing  velocity  v<\ep  (m  s-1)  directed  in  the  normal  di¬ 
rection  of  the  electrode  into  the  electrode  domain,  can  be  deduced 
from  the  molar  flux,  the  density  and  the  molar  mass  of  zinc. 


^dep 


MznNzn 

Pin 


(12) 


2.2.3.  Boundary  conditions 

In  this  work,  zinc  electrodeposition  model  was  constructed  us¬ 
ing  COMSOL  Multiphysics  software  and  kinetics  Monte  Carlo 
method.  To  simplify  the  model,  we  assume  that  the  electrolyte  flow 
with  laminar  and  similarity  technique  is  applied  to  surface  rough¬ 
ness  of  the  cathode.  Different  flowing  velocities,  different  charging 
modes,  and  different  surface  morphology  of  the  electrode  were 
investigated  on  the  basis  of  the  model. 

A  schematically  three-dimensional  model  for  the  zinc  electro¬ 
deposition  is  shown  in  Fig.  2,  which  contains  four  governing 


3.1.  Zinc  regeneration  process 

Zinc  electrodeposited  process  includes  three  steps,  namely  ionic 
transport,  electrochemical  reaction  and  electrocrystallization, 
involving  mass  transfer,  charge  transfer,  and  zinc  adatoms  incor¬ 
porated  into  metal  lattice.  The  depositing  morphology  at  the  sur¬ 
face  of  the  cathode  is  mainly  controlled  by  reaction  kinetics, 
diffusion,  atomic  bond,  crystal  nucleation,  and  crystal  growth  in 
space  and  time.  The  morphology  of  electrodeposited  zinc  is  the 
balance  result  of  thermodynamics  and  kinetics,  as  illustrated  in 
Fig.  3,  thereby  bringing  about  three  different  growth  modes 
including  the  Volmer-Welmer  (VW)  mode,  the  Frankvan-der- 
Merwe  (FM)  mode,  and  the  Stranski-Krastanov  (SK)  mode  [23]. 
Zinc  atoms  uniformly  grow  at  the  horizontal  direction  controlled  by 
diffusion  along  the  depositing  surface,  while  mass  transfer  between 
atomic  layers  controlled  by  Eheilich-Schwoebel  barrier  (Es) 
determine  the  uniformity  at  the  vertical  direction.  The  final 
morphology  depends  on  the  competition  between  the  vertical 
growth  and  the  horizontal  growth  of  crystal  planes. 

Non-uniform  morphology  of  deposited  zinc  may  be  induced  by 
crystal  imperfections  such  as  point  defects,  edge  and  screw  dislo¬ 
cations  and  crystal  clusters  [7  .  Fig.  4  shows  several  depositing 
types  of  adatoms  and  islands  on  the  substrate  at  the  initial  stage 
and  the  shape  change  of  deposited  zinc  at  the  later  stage  from  the 
atomic  perspective.  The  morphological  evolution  of  electro¬ 
deposited  zinc  may  be  simulated  by  means  of  kinetics  Monte  Carlo 
method  [24  ,  where  zinc  atom  number  at  each  depositing  layer  is 
unequal  and  decreases  progressively  due  to  energy  disparity  of  zinc 
atoms  by  van  der  Waals'  force  and  different  overpotentials  at 
different  positions,  resulting  in  morphological  change. 

The  morphology  of  deposited  zinc  can  be  in  the  form  of  spongy, 
dendritic,  boulder-like,  layer-like,  mossy,  and  filamentous  [7], 
mostly  depending  on  two-phase  interface  conditions  and  external 
charging  mode.  Zinc  dendritic  growth  controlled  by  diffusion  is  one 
of  morphological  change,  which  can  be  described  by  the  method  of 
diffusion-limited  aggregation  [25].  Accordingly,  many  researchers 
described  dendritic  morphology  and  characteristics  in  metal  elec¬ 
trodeposition  on  basis  of  fractal  theory  26,27  .  Layer-like  structure 
arises  from  epitaxial  growth  along  with  crystal  structure  and 
orientation.  Boulder-like  deposit  is  known  as  isoprobability  growth 
at  different  crystal  faces  and  orientations.  Mossy  deposit  is 
controlled  by  activation.  These  morphologies  are  closely  linked  to 
electrode  structure,  electrolyte  hydrodynamics,  and  charging 
mode. 

3.2.  Structure  optimization 

Zinc  electrode  is  essential  in  zinc  regeneration  either  for 
zinc-air  batteries  or  zinc-air  fuel  cells  [28  ,  where  zinc  electrode 
structure  as  a  substrate  of  deposited  zinc  influences  morphological 
evolution. 

Zinc  depositing  surface  was  assumed  to  be  comparatively  uni¬ 
form,  not  considering  surface  roughness  of  zinc  electrode  such  as 
the  planar  electrode.  Morphological  changes  of  deposited  zinc 
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Table  2 

Variables  and  equations. 


Quantity 

Variable 

Equation 

Position 

(x,  y,  z ) 

L  =  a  (m),  W  =  b  (m),  H  =  d  (m) 

Electrolyte  velocity 

V 

v  =  QJbd  (m  s”1) 

Peclet  number 

Pe 

Pe  =  dv/D 

Flux  density 

(p 

q>  =  RT/F  (V) 

occur  at  the  zinc  electrode  surface  during  charge,  as  shown  in 
Fig.  5(a),  where  more  zinc  is  deposited  at  the  cathode  edges  than  at 
the  cathode  center,  but  an  alternative  presents  the  preferable  result 
by  means  of  downsizing  the  cathode  and  smoothing  the  edges  of 
the  cathode  (Fig.  5(b)).  This  phenomenon  may  be  attributed  to  edge 
effect  of  the  zinc  electrode,  leading  to  polarization  easily  at  the 
periphery  of  the  electrode.  Therefore,  the  cathode  and  the  anode 
should  be  kept  parallel  in  structure,  and  the  dimension  of  the  anode 
is  larger  than  that  of  the  cathode,  among  which  the  cathode  edges 
need  to  be  round  off. 

Increasing  surface  roughness  of  zinc  electrode  like  a  porous 
electrode  will  cause  morphological  change,  where  the  morphology 
of  deposited  zinc  would  be  uneven  over  the  surface  of  electrode  at 
the  early  depositing  stage,  namely  the  more  zinc  deposited  at  the 
bulges,  the  less  at  the  scallops.  Fig.  6  shows  local  current  densities 
distributed  on  the  rough  surface  of  the  cathode,  where  the  local 
current  densities  is  higher  at  the  bulges  than  at  the  scallops. 
Consequently,  different  surface  microstructures  lead  to  different 
current  densities  at  the  different  points  of  the  electrode,  and  then 


Fig.  3.  Growth  mode  of  electrodeposited  zinc. 


induce  different  deposit  velocities.  Besides,  the  porous  electrode 
with  high  active  area  is  available  for  diffusion  of  depositing  atoms. 

When  an  electrode  of  discrete  columnar  structure  is  employed 
as  the  cathode,  the  morphology  of  deposited  zinc  tends  to  be 
granular  analogous  to  the  column  dimension.  Fig.  7  shows  the 
morphology  of  deposited  zinc  at  the  direct  current  density  of 
100  mA  cm-2,  where  the  depositing  mass  change  at  the  edges  are 
higher  than  at  the  centers  due  to  edge  effects  of  electric  field.  It  is 
even  more  important  that  zinc  atoms  only  deposit,  diffuse  and 
aggregate  on  the  small  electrode  surface,  thus  strengthening 
cohesion  force  among  zinc  atoms. 

The  gas  bubbles  produced  by  electrochemical  reaction  play  a 
predominant  role  in  the  depositing  morphology  as  well.  Zinc-air 
battery  was  placed  as  the  vertical  position,  leading  to  aggregation 
of  oxygen  bubbles  at  the  upper  electrolyte  and  driving  the  depos¬ 
iting  atoms  flowing  towards  downside  during  charge  process,  as 
shown  in  Fig.  8(a).  Zinc-air  battery  was  placed  as  the  horizontal 
position,  timely  removing  oxygen  bubbles  produced  by  reduction 
reaction,  and  the  compact  and  uniform  deposit  layer  was  attained 
in  Fig.  8(b). 

3.3.  Electrolyte  management 

Zincate  concentration  in  the  electrolyte  affects  the  depositing 
morphology  as  well.  Zincate  concentration  gradient  will  be  formed 
between  the  bulk  electrolyte  and  the  interfacial  electrolyte  near  the 
cathode  as  a  result  of  electrochemical  reaction,  and  zincate  con¬ 
centrations  at  different  reaction  zones  will  be  different  due  to 
different  reaction  rates.  For  the  purpose  of  uniformly  stable  species 
concentration  of  solution,  it  is  necessary  to  make  the  electrolyte 
flow. 

Assuming  that  the  electrolyte  move  by  way  of  laminar  flow,  the 
electrolyte  between  two  parallel  electrodes  flows  from  one  side 
towards  the  other  side,  with  the  middle  part  of  the  flowing  velocity 
is  higher  than  top  and  bottom  parts.  Now  that  convection  induced 
by  electrolyte  hydrodynamics  would  hold  dominant  in  species 
transport  in  the  electrolyte  relative  to  diffusion,  there  is  a  need  to 
take  full  account  of  ions  transport,  electrolyte  concentration  and 


Fig.  4.  Schematic  process  of  electrodeposited  zinc. 
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Surface:  Deposited  zinc  mass  distribution^  rrf:) 
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Fig.  5.  Mass  distribution  of  deposited  zinc  of  parallel  electrodes  with  different  structures. 
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Surface:  Deposited  zinc  mass  distribution^  m‘:) 
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dendritic  growth  to  optimize  the  flowing  field.  The  flowing  velocity 
of  electrolyte  can  affect  the  distribution  of  species  concentration 
and  the  morphology  of  zinc  deposition.  The  larger  the  velocity,  the 
more  uniform  the  concentration,  as  shown  in  Fig.  9.  However,  the 
high  velocity  will  scour  nearby  layer  of  deposited  zinc  owing  to  the 
small  binding  force  among  atoms  on  the  sediment  surface,  and 
then  cause  shape  change. 

Fig.  10  shows  that  two  different  zinc  depositing  morphologies 
are  obtained  under  the  conditions  of  the  quiescent  electrolyte  and 
the  flowing  electrolyte.  Dendrites  did  not  appear  mainly  because  of 
convective  flow  of  the  electrolyte  reducing  concentration  gradients 


between  bulk  electrolyte  and  surface  electrolyte  near  the  cathode, 
and  taking  away  depositing  atoms  around  dendrites.  It  is  proved 
that  the  flowing  electrolyte  is  an  effective  method  to  solve  the 
problem  of  dendritic  growth.  However,  the  flowing  electrolyte  also 
washes  away  zinc  deposits  near  the  ingress  and  aggregated  on 
other  sites. 

The  flowing  electrolyte  can  promptly  take  away  gas  bubbles 
avoiding  gas  aggregation  near  the  anode.  However,  the  high  flow¬ 
ing  velocity  will  adversely  affect  species  transport,  thus  affect  local 
current  density  of  zinc  electrode.  Fig.  11  shows  the  local  current 
density  for  the  four  Peclet  numbers.  The  higher  flow  velocity  for  the 


Surface:  Local  current  density  (Am^) 
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Fig.  6.  Local  current  densities  on  the  coarse  surface  of  zinc  electrode. 
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Surface:  Electrode  mass  change  (kg  m'2) 
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Fig.  7.  Deposited  zinc  morphology  on  the  discrete  pins  of  electrode. 


higher  Peclet  number  will  increase  the  local  current  density 
significantly  due  to  the  increased  transport  velocity  of  zinc  ions. 

It  is  found  in  Fig.  12  that  maximum  thickness  of  deposited  zinc 
decrease  under  the  condition  of  different  fluid  velocities,  where 
increased  flow  rate  cannot  limit  dendrites  growth  at  small  flowing 
velocities,  properly  because  of  this  flow  rate  ameliorating  concen¬ 
tration  polarization  between  bulk  electrolyte  and  reaction  zone  and 
improving  local  current  densities  on  the  cathode  surface.  However, 
zinc  deposited  thickness  decreases  with  increased  flow  rate, 
attributed  to  driving  force  of  hydrodynamic  electrolyte  exceeding 
bonding  energy  among  zinc  atoms. 

In  view  of  the  above  mentioned,  it  is  found  that  the  flowing 
electrolyte  not  only  mitigates  the  instability  of  electrolyte  charac¬ 
teristics  during  electrodeposition,  but  electrolyte  hydrodynamics 
helps  to  improve  ion  diffusion  and  suppress  zinc  dendritic  growth 


and  take  away  gas  bubbles  in  the  zinc-air  battery  with  alkaline 
electrolyte. 

3.4.  Charging  mode 

It  is  well-known  that  an  increase  in  charging  current  accelerates 
the  rate  of  ion  migration,  increases  the  number  of  activated  ions, 
and  enhances  atom  energy.  Activation-controlled  deposition  plays 
a  leading  role  at  low  current  densities  due  to  ions  translating  into 
atoms  overcoming  the  activation  energy  barrier,  while  diffusion- 
controlled  deposition  dominates  at  high  current  densities.  Fig.  13 
shows  that  depositing  morphology  is  the  same  as  boulder-like  at 
high  current  density,  while  compact  at  low  current  density. 

Fig.  14  shows  the  trend  of  maximum  deposited  zinc  thickness  at 
different  charging  current  densities  of  100  A  nrT2-2000  A  m-2  at 


Fig.  8.  Deposited  zinc  morphology  in  weight  40%  KOH  with  saturated  zinc  oxide  for  10  min  at  the  current  density  of  100  mA  cm  2:  (a)  vertical  position,  (b)  horizontal  position. 
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Pe  =  100  Surface:  Concentration  (mol  m'2) 
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Fig.  9.  Electrolyte  concentration:  (a)  Pe  =  100,  and  (b)  Pe  =  le5. 

the  electric  quantity  of  3600  A  s.  The  thickness  of  deposited  zinc 
displays  as  a  positive  relationship  with  charging  current  mainly  due 
to  more  electric  quantity  for  zinc  ions,  therein  the  overpotential  of 
zinc  oxide  reduction  reaction  increases  towards  minus  direction 
resulting  from  polarization,  and  the  voltage  between  the  anode  and 
the  cathode  increases  with  current  densities.  Zinc  atoms  will 
quickly  diffuse  and  grow  on  the  electrode  surface  at  high  current 
densities,  leading  to  shorting  and  failure  of  zinc-air  batteries. 


Meanwhile  much  more  electrons  may  be  wasted  for  hydrogen 
formation  or  non-Faradic  reaction.  This  means  that  zinc  regenera¬ 
tion  in  zinc-air  battery  need  to  take  long  time  with  small  charging 
current,  thus  zinc  morphological  control  at  high  charging  currents 
is  crucial  for  zinc-air  batteries  in  application  of  energy  storage. 

Phase  changes  occur  next  to  the  electrodes  during  electrode¬ 
position,  including  liquid-solid  phase  change  and  liquid-gas  phase 
change.  The  production  of  hydrogen  is  companied  by  zinc  elec¬ 
trodeposition  which  depends  on  electrolyte  and  electrode 
composition,  charge  mode,  and  temperature.  Hydrogen  bubbles 
can  give  birth  to  electrolyte  convection,  and  thus  affect  zinc 
redistribution  over  the  electrode.  It  is  worth  noting  that  gas  bubbles 
may  stay  within  depositing  layer,  and  then  lead  to  shape  change 
especially  for  planar  electrode,  as  shown  in  Fig.  15.  To  reduce 
hydrogen  evolution,  it  is  required  to  add  a  few  inorganic  elements 
or  organic  compounds  into  the  electrolyte  to  increase  the  hydrogen 
evolution  overpotential  [29,30  . 

In  addition,  growth  of  dendrites  is  mainly  controlled  by  diffu¬ 
sion  at  high  current  densities.  Pulsating  current  or  periodic  reverse 
current  deposition  is  conducive  to  improve  the  depositing 
morphology  due  to  zincate  ions  having  enough  diffusion  time. 
Fig.  16  displays  that  distinct  results  in  morphology  and  depositing 
mass  were  obtained  under  conditions  of  direct  current  and  pul¬ 
sating  current,  respectively.  It  is  found  that  the  depositing 
morphology  of  zinc  is  greatly  affected  by  the  charging  mode  with 
duty  cycle  which  is  defined  as  Tonl(Ton  +  T0ff).  Reducing  the  duty 
cycle  effectively  improve  the  morphology  of  deposited  zinc,  namely 
compact  morphology  at  the  macrostructure  scale  and  layer-like 
morphology  at  the  micrometer  scale  attained  under  the  condition 
of  pulsating  currents.  Therefore,  appropriate  charging  mode  is  an 
effective  method  to  suppressing  the  dendritic  growth. 

Pulsating  currents  can  improve  the  morphology  of  deposited 
zinc,  namely  pulse  current  with  different  duty  cycles.  Maximum 


Fig.  10.  Deposited  zinc  morphologies  in  weight  40%  KOH  with  1  M  zinc  oxide  with  electric  quantity  for  10  min  at  the  current  density  of  100  mA  cm  2,  (a)  quiescent  electrolyte,  (b) 
flowing  electrolyte. 
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Fig.  11.  Local  current  densities  on  the  electrode  surface  for  different  Peclet  numbers. 
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Fig.  12.  Morphological  evolution  of  deposited  zinc  at  different  flowing  velocities  of  the  electrolyte. 


Fig.  13.  Zinc  depositing  morphology  in  weight  40%  KOH  with  saturated  zinc  oxide  for  10  min  under  vertically  positioned  condition:  (a)  50  mA  cm  2  of  current  density,  (b) 
20  mA  cm  2  of  current  density. 
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Fig.  14.  Morphological  evolution  of  deposited  zinc  at  different  charging  currents. 


thickness  of  deposited  zinc  tends  to  decrease  with  duty  cycle 
reduction  at  the  electric  quantity  of  16.2  A  s  in  Fig.  17,  as  a  result  of 
diffusion  of  zinc  atoms  for  off-time  cycle.  Therefore,  morphological 
evolution  is  under  control  by  means  of  charging  mode,  where  high 
current  combined  with  pulsating  mode  or  small  direct  current  is 
available  for  uniform  morphology  of  deposited  zinc. 

Based  upon  the  above  results,  we  can  see  that  shape  change  of 
zinc  morphology  on  the  cathode  is  closely  related  to  local  current 
densities  under  the  same  external  conditions.  However,  dendritic 
growth  may  be  effectively  controlled  by  variables  such  as  flowing 
electrolyte  or  pulsating  current.  Additionally,  oxygen  or  hydrogen 
bubbles  not  only  directly  influence  the  depositing  morphology,  and 
deteriorate  zinc-air  battery's  performance. 

4.  Conclusion 

The  mechanism  for  the  morphological  change  of  zinc  electro¬ 
deposition  in  a  zinc-air  battery  was  investigated  by  means  of 


three-dimensional  kinetic  Monte  Carlo  method  and  COMSOL  soft¬ 
ware.  The  morphologies  of  deposited  zinc  under  various  testing 
conditions  were  also  studied  by  way  of  scanning  electron  micro¬ 
scopy.  The  results  show  that  the  morphological  change  of  deposited 
zinc  mainly  depends  on  local  current  densities  of  the  cathode  and 
are  closely  linked  to  working  conditions  as  well,  including  the  hy¬ 
drodynamic  electrolyte,  charging  mode,  surface  quality  of  the 
electrode,  and  structure.  The  detail  is  as  follows: 

1)  Zinc  morphology  development  depends  on  the  rates  of  ion 
transfer,  charge  transfer  and  atom  attachment  at  different 
depositing  conditions,  stemming  from  the  factors  of  diffusion, 
activation  and  electrocrystallization.  Different  deposited  mor¬ 
phologies  are  caused  by  different  determining  factors,  among 
which  dendritic  growth  is  mainly  controlled  by  diffusion. 

2)  Zinc  deposited  morphology  control,  where  low  current  den¬ 
sities,  duty  cycle  of  pulsating  currents  and  flowing  electrolyte 
are  in  favor  of  mass  and  charge  diffusion  as  well  as  high  zincate 


Fig.  15.  (a)  Zinc  depositing  morphology  at  50  mA  cm  2  of  current  density,  (b)  distribution  of  local  current  densities. 
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Fig.  16.  Zinc  depositing  morphology  with  electric  quantity  for  10  min:  (a)  100  mA  cm  2 


of  direct  current  density,  (b)  100  mA  cm  2 


of  pulsating  current  density,  20%  duty  cycle. 
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Fig.  17.  Morphological  evolution  of  deposited  zinc  at  different  duty  cycles. 
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concentration  and  low  KOH  concentration,  thus  resulting  in 
compact  morphology  of  deposited  zinc. 

3)  Even  distribution  of  deposited  zinc  on  the  cathode  surface, 
where  horizontal  layout  of  electrodeposited  device  is  a 
preferred  design  for  evenly  distributed  zinc  deposition  relative 
to  vertical  layout.  In  addition,  electrolyte  hydrodynamics  and 
edge  effect  of  zinc  electrode  also  cause  uneven  distribution  of 
deposited  zinc. 

4)  Phase  changes  of  zinc,  hydrogen  and  oxygen  during  charge, 
where  these  phase  changes  influence  morphological  evolution 
and  distribution.  Hydrogen  evolution  can  cause  convective 
electrolyte  to  increase  charge  transfer  rate,  resulting  in 
randomly  nucleation,  and  oxygen  aggregation  will  separate  the 
anode  from  the  electrolyte,  leading  to  reduction  of  current 
efficiency. 
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